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Abstract — It is well known that opportunistic scheduling algo- 
rithms are throughput optimal under full knowledge of channel 
and network conditions. However, these algorithms achieve a 
hypothetical achievable rate region which does not take into ac- 
count the overhead associated with channel probing and feedback 
required to obtain the full channel state information at every slot. 
We adopt a channel probing model where /3 fraction of time slot 
is consumed for acquiring the channel state information (CSI) 
of a single channel. In this work, we design a joint scheduling 
and channel probing algorithm named SDF by considering the 
overhead of obtaining the channel state information. We first 
analytically prove SDF algorithm can support 1 + e fraction of 
of the full rate region achieved when all users are probed where e 
depends on the expected number of users which are not probed. 
Then, for homogenous channel, we show that when the number 
of users in the network is greater than 3, e > 0, i.e., we guarantee 
to expand the rate region. In addition, for heterogenous channels, 
we prove the conditions under which SDF guarantees to increase 
the rate region. We also demonstrate numerically in a realistic 
simulation setting that this rate region can be achieved by probing 
only less than 50% of all channels in a CDMA based cellular 
network utilizing high data rate protocol under normal channel 
conditions. 

I. Introduction 

Scheduling is an essential problem for any shared resource. 
The problem becomes more challenging in a dynamic setting 
such as wireless networks where the channel capacity is time 
varying due to multiple superimposed random effects such as 
mobility and multipath fading. Opportunistic scheduling has 
emerged as an attractive solution for improving the efficiency 
of wireless systems with limited resources such as frequency 
band and power. In principle, opportunistic policies schedule 
the user with the favorable channel conditions to increase the 
overall performance of the system fl/]. Optimal scheduling in 
wireless networks has been extensively studied in the literature 
under various assumptions. The seminal work by Tassiulas 
and Ephremides have shown that a simple Lyapunov-based 
opportunistic algorithm that schedules the user with the highest 
queue backlog and transmission rate product at every time slot, 
can stabilize the network, whenever this is possible J2). 

There has been much work in developing scheduling algo- 
rithms for down-link wireless systems for various performance 
metrics that include stability, utility maximization and energy 
minimization Q, HI, Q. However, the common assumption 
in these studies is that the exact and complete channel state 
information, (CSI) of all users is available at every time slot. 
Hence, these algorithms achieve a hypothetical rate region 
by assuming that full channel state information is available 
without any channel probing or feedback costs. However, 
in practice acquiring CSI introduces significant overhead to 



the network, since CSI is obtained either by probing the 
channel or via feedback from the users. In current wireless 
communication standards such as WiMax [4] and LTE 0, 
there is a feedback channel used to relay CSI from the users to 
base station. Obviously, this feedback channel is bandlimited 
and it is impossible to obtain CSI from all users at the same 
slot. As a motivating example, consider CDMA/HDR (High 
Data Rate) system J6|, where the Signal-to-Noise Ratio (SNR) 
of each link is measured. The value of the SNR is then 
mapped to a value representing the maximum data rate that 
can support a given level of error performance. This channel 
state information is then sent back to the base station via the 
reverse link data rate request channel (DRC). The channel state 
information is updated every 1.67 ms and sent back as 4 bits. 
Assume that there are 25 users in a cell, and hence, 100 bits of 
channel information has to be sent every 1.67 ms. This requires 
60Kbps of channel rate to be dedicated only for channel 
measurements. The minimum data rate of HDR system is 38.4 
Kbps and the average data rate is 308Kbps. Thus, overhead 
of acquiring CSI is twice the minimum data rate, and is 
approximately more than 20% of the average transmission 
rate. This overhead increases significantly in a multichannel 
communication system such as LTE. Thus, channel probing 
must be done efficiently in order to balance the trade-off 
between being opportunistic (e.g., obtaining useful channel 
information) and consuming valuable resources. 

Most of real systems run under limited hardware and 
software capacity. Hence, in addition to taking into account 
the probing cost, the implementability of feedback schemes 
must be ensured in practical systems. For instance, a practical 
feedback scheme must not require high computation time 
and overhead for probing users. In addition, any feedback 
algorithm must not use any statistics such as channel or 
arrival statistics for probing decision since in real word, such 
statistics change over time. It is also beneficial to design a 
limited feedback system that can operate over a wide channel 
condition and numbers of users. 

In this work, we consider in fully connected network (e.g., 
Cellular network, WLAN) where base station (BS) is trans- 
mitting to a fixed number of users. We assume that each user 
has infinite buffer capacity and data arrives into users' queues 
according a stochastic arrival process. We aim to develop a 
joint scheduling and channel probing algorithm that stabilizes 
users' queues with taking into account to the probing cost. Our 
algorithm dynamically determines the set of channels that must 
be probed at every time slot based on the information obtained 
from the channel state information of the user that has the 



maximum queue length at a given time slot. The scheduling 
part is based on well known Max- Weight algorithm J2]. 
Our contributions are summarized as follows: 

• We first propose a joint scheduling and channel probing 
algorithm which is easy to implement and low cost. In 
addition, our algorithm does not require any statistics 
(i.e., channel or arrival statistics) and can be applied over 
correlated and even non-stationary channels. 

> Technical contributions of our paper are as follows: we 
first show that the proposed algorithm can achieve 1 + e 
fraction of the full rate region of the case when all users 
are probed. For homogenous channels, we analytically 
show that e > (e.g., we guarantee to increase the rate 
region) when the number of users is greater than 3. For 
heterogenous channels, we prove that e > as the number 
of channel states goes infinity or the number of users is 
large enough. 

• We implement a realistic network setting where we sim- 
ulate High Data Rate (HDR) protocol in CDMA cellular 
networks and show by numerical analysis that when our 
proposed algorithm is used a comparable performance 
with Max- Weight algorithm with full CSI can be achieved 
by only probing on the average less than 50% of users. 

The rest of the paper is organized as follows. Section HH 
summarizes the literature on opportunistic scheduling algo- 
rithms considering probing overhead. Section |TTT] presents the 
network model, the basic structure of the channel probing 
model. In Section IIVI we give the problem formulation under 
the proposed channel probing models. In Section [IV] we give 
joint scheduling channel probing algorithms for both channel 
probing models. Numerical results are presented in Section 
IVIII Finally, Section IVIIII concludes the paper and presents 
possible future directions. 

II. Related Works 

There has been a significant recent interest in applying 
Max- Weight-like algorithms in more realistic wireless network 
settings. Throughput-maximizing scheduling has been studied 
with different forms of limited CSI. For instance, infrequent 
channel state measurements was investigated in [7] and it was 
shown that achievable rate region shrinks as the frequent of 
CSI decreases. The impact of delayed CSI was investigated 
in 0. However, unlike these works, we aim to minimize the 
number of probed users at a time while maximizing the achiev- 
able rate region. In this context, one of the most intriguing 
research challenges in the context of wireless networking with 
limited feedback is the design of a scheduling policy that 

1) is implementable, simple and low-complexity, 

2) achieves high performance, i.e., low packet delay and 
large rate region. 

3) works (correlated or even over non-stationary channels) 
without requiring any statistics such as channel distribu- 
tion. 

In the following, we classify the works along this direction. 
Limited feedback bandwidth: In (9), the authors proposed 



a joint scheduling and channel probing algorithm stabilizing 
the network by allowing the base station to probe a subset of 
channels (or links) at each time slot. First, the a throughput- 
optimality algorithm was developed that uses the expected 
channel rates, i.e., the channel distributions are given. While 
it is reasonable to estimate the joint channel state distribution 
when channels are independent, when the number of channels 
is small and channel statistics do not change. However, such 
estimation becomes intractable in real word where channels 
are non-stationary process and when the number of channels 
is large. Then, the authors proposed only queue length based 
algorithm. However, throughput-optimality of that algorithm 
can only be shown under certain condition, i.e., when chan- 
nels are symmetric and subsets of channels are disjoint. In 
addition, delay performance of these algorithms are unknown. 
In iflOl . the authors proposed a variant of the algorithm in [9] 
to analyzes queue-overflow performance for scheduling with 
limited CSI. 

In ifTTI . a feedback allocation algorithm was proposed for 
multi-channel system with limited feedback bandwidth. In 
other words, only a limited number of users can be probed at 
a time. It was shown that the proposed algorithm can achieve 
1 — e fraction of the full rate region when channel distributions 
are known at the BS. Note that the algorithms in J9], (TT], [ED 
cannot achieve all three properties given above. 

Unlike the works in |9), IfTTI . IflOl . we assume a more 
flexible channel probing model where BS may probe all 
users at a given time. In other words, we did not assume 
bandwidth limited feedback channel. Similar model was used 
in iTPH . However, it was assumed that all users have saturated 
buffers, i.e., they always have packets in the corresponding 
buffers and, hence, the network stability problem was not 
investigated. In this sense, the most related work is iTPJI . 
where the optimal feedback- scheduling scheme for a 
single-channel downlink is derived. Specifically, in iTPJl . it 
was assumed that a single channel probing requires a certain 
portion of data transmission (i.e., f3 fraction of data slot). 
Hence, the server has a cost for probing channels and gains 
a reward (queue weighted throughput) which depends on 
the user. The problem of finding optimal joint algorithm is 
transformed into an optimal stopping time problem and is 
solved by Markov Decision Process (MDP) where the channel 
probabilities are known to the BS. In addition, since the 
authors uses MDP to solve the problem, this formulation is 
computationally intractable as it involves a high dimensional 
state. Therefore, this work does not satisfy all three properties 
either. 

Markovian Channels: Studies in fl4l. ifBl. lfl6l. ATI. ifTsl 
attempt to learn the channel distribution for scheduling. 
In lfl8ll . a 2-stage decision procedure is used where channel 
measurement and packet transmission has to be performed 
at every time slot without knowing the channel distribution. 
This work attempts to learn transmission success probabilities 
by averaging the previous outcomes. Similarly, in lfl4l . the 
authors proposed to estimates the channel statistics by using 
some portion of the time slots for observation slot with 



some probability. In |[T5l , the authors proposed a scheduling 
algorithm under imperfect CSI in single-hop networks with 
i.i.d. channels. They consider that probing a channel brings 
a certain amount of energy cost. Under this setup, the 
proposed scheduling algorithm which decides whether to 
probe the channel with the energy cost is a Max- Weight type 
scheduling policy that minimizes the energy consumption 
subject to queue stability. In ifTTll . it was assumed that 
wireless channels evolve as Markov-modulated ON/OFF 
processes. With this assumption, a exploitation-exploration 
trade-off was investigated. Similarly,in [16], a two state 
discrete time Markov chain with a bad state which yields no 
reward and a good state which yields reward was considered. 
Aforementioned works assume that the underlaying stochastic 
process of the channel evolves according to a fixed stationary 
process such as ergodic Markov chain. In practice, such an 
assumption does not hold most of the time. For instance, the 
measurement study in [19| shows that the wireless channel 
exhibits time-correlated and non-stationary behavior. Thus, 
these works cannot achieve all three properties given above. 

III. System Model 

We consider a cellular system with a single base station 
transmitting to N users with a fixed power. Let N denote 
the set of users in the cell. Time is slotted, t G {0, 1,2,.. .}, 
and wireless channel between the base station and a mobile 
user is assumed to be independent across users and slots. 
The gain of the channel is constant over the duration of a 
time slot but varies between slots. In practical systems (e.g., 
CDMA/HDR system 0), transmission rate is determined by a 
link adaptation algorithm, which selects the highest transmis- 
sion rate to meet a given allowable target error probability. 
Only finite number of transmission rates can be supported 
due to modulation and coding schemes. We assume that 
each channel has L possible states with corresponding rates 
1Z = {n, r%, ■ ■ ■ , tl] listed in descending order, i.e., r^ > n 
if I > k. Note that r^, k G {1,2, ...,L} only depends on 
the Signal-to-Noise Ratio (SNR). We denote R n (t) G 1Z as 
the channel state information (CSI) of user n at time t. For 
user n, let p£ denote the probability of being state k, i.e, 
Pr[R n (t) = r*.] = Pfr and each channel state is possible with 
non-zero probability such that p mm < p™ < p max , \/ n ,i. 

A. Channel Probing Model 

In this work, we assume that there is no dedicated feedback 
channel, and CSI is relayed over the data channel. Hence, 
depending on the scheduling algorithms, the scheduler can 
probe any number of users at any time slot. We quantify the 
overhead of obtaining the CSI of a single user in terms of a 
time fraction of the time slot as in lfl31 . This time duration 
may include the time spent for pilot signal transmission, 
measurement of the signal strength of pilot signal and the 
transmission of CSI to the base station. 

We assume that at the beginning of a time slot, the BS sends 
a pilot signal. Based on the quality of the received pilot signal, 



each user n can determine its current channel state R n (t). We 
also assume that downlink and uplink channels are identical, 
i.e., the maximum transmission rates in both directions are the 
same. Hence, the BS obtains reports from a selected number of 
users before scheduling a downlink transmission. Let N p (t) C 
M be the set of users who report back their CSI. Let I n (t) be 
an indicator function showing whether user n is scheduled to 
be transmitted to in slot t: 



In(t) 



1 ; if user n G Af p (t) and it is scheduled 
; otherwise 



(1) 



Note that for successful transmission, a user scheduled to be 
transmitted to in slot t should be selected from among the 
users who have reported their channel states. 

We assume that j3 fraction of the time slot is consumed to 
obtain CSI from a single user. Hence, only (1 — f3N p (t)) x T s 
seconds are available for data transmission where N p {t) is the 
cardinality of N p (t) and T s is the duration of the time slot. 
We assume that 1 - PN p (t) < INote that when full CSI is 
obtained, this time is equal to (1 — f3N) x T s , We assume that 
at most one user can be scheduled at a given time. Hence, the 
amount of data that can be transmitted to user n by the BS 
when N p (i) users are probed at time t is given by, 



D n (t) = (1 - (3N P (t))T s R n {t)I n {t). 



(2) 



We assume that T s is normalized to unit slot length, i.e., T s = 
1 in the rest of the paper. Let A n (i) be the amount of data 
(bits or packets) arriving into the queue of user n at time 
slot t. We assume that A n (t) is a stationary process and it is 
independent across users and time slots. We denote the arrival 
rate vector as A = (Ai, A2, . • . , Ajv), where A„ = E[A„(i)]. 
Let Q(t) = (Qi(t),Q 2 (t), . . . ,Q N (t)) denote the vector of 
queue sizes, where Q n (t) is the queue length of user n at 
time slot t. 

Definition 1: A queue is strongly stable if 



1 



t—x 



limsupy y ]E[Q n (t)} < 00 

T = 



(3) 



Moreover, if every queue in the network is stable then the 
network is called stable. The dynamics of the queue of user 
n is given by, 

[Q n (t) + A n (t) - D n (t)} + . 



»(<+!) 



(4) 



where [x} + = max(x, 0). 

Before giving the proposed scheduling algorithm the fol- 
lowing definitions are useful. 

The achievable rate region or rate region of a network is 
defined as the closure of the set of all arrival rate vectors A 
for which there exists an appropriate scheduling policy that 
stabilizes the network. 

Definition 2: Ah is the hypothetical rate region where full 
CSI is available (e.g. by an Oracle) without any channel 
probing or feedback costs, i.e., j3 = 0. 

As discussed in ll20l . it is impossible to achieve the boundary 
of Ah in real systems, since there is an overhead of acquiring 
CSI from the users. 



Definition 3: Af is the achievable rate region when probing 
cost is taken into account and when all users' channels are 
probed at every time slot according to the feedback model. 
Note that A/ c A h 

Definition 4: A a is the algorithm-dependent achievable rate 
region where the probing cost is taken into account. 

Note that if Af C A a C A/,, then we can argue that there 
exists a solution or algorithm which is more efficient than full 
CSI Max- Weight algorithm in terms of the rate region. We 
define the weighted rate of user n as follows: 



W n (t) = Q n (t)R n (t), 



(5) 



and w(t) is the maximum weighted rate such that w(t) = 
argmax rlGA/ - {W„(4)}. In the next section, we propose a 
scheduling and dynamic feedback algorithm which finds the 
maximum weighted rate at every time slot by probing only 
N p (t) number of users where N p (i) < N. 

IV. Scheduling and dynamic Feedback (SDF) 
Algorithm 

Recall that our aim is to find a joint scheduling and channel 
probing algorithm which is able to find the user which has 
the maximum weighted throughput at each time slot with 
minimum number of channel probing. In this section, we 
propose a joint scheduling and dynamic feedback allocation 
algorithm that determines the user which has the maximum 
weighted throughput at each time slot without probing every 
user. 



SDF Algorithm: 

(1) probing decision: 

• Step 1: Determine the user which has the maximum 
queue length, 

i* = argmax{<5j(i)} 

• Step 2: Probe and acquire the CSI of user i*. Let 
Ri* (t) be the CSI of user i* at time t. 

• Step 3: Broadcast the value of Ri*(t). 

• Step 4: The users which have higher rate than Ri* (t) 
report their CSIs to BS. 

Let us define the set S p (t) as follows: 

S p (t)±{j:R j (t)>R i ,(t)} 

(2) scheduling decision: 

BS schedules a user which has the maximum weighted 
throughput according to Max- Weight algorithm J3] as 
follows; 

n* = argmax{(l - f3N p (t))W n (t)} (6) 



where 



Af p (t)£-S p (t)Ui* 



Intuition: Let us assume that user n* has the maximum 
weighted rate at a time slot. Given CSI of the user i*, the 
CSI of users with rate lower than Ri* need not be collected 
since their weighted rates are always smaller than that of user 
i*. SDF algorithm is especially efficient when the number of 
users is large since in that case the number of users with rate 
lower than that of user i* is large with high probability. 

A. Analysis of SDF Algorithm with Homogenous Channels 

Now, we investigate the increase in achievable rate region 
when SDF algorithm is employed. We first consider a homoge- 
nous channel model where, 

Pk =Pk,Vn. 

Let M(t) denote the number of users which do not send 
their CSIs since their channel conditions are worse than the 
user which has the maximum queue length at time t. The 
number of fractions of time slot consumed for probing with 
SDF algorithm is determined as follows: first, the BS acquires 
the CSI of user i* and /3 fraction of time slot is used for 
probing user i* . Then, BS broadcasts the value of CSI of user 
i*. We assume that broadcasting CSI of a user also consumes 
/3 fraction of time slot. Then, the number of users which have 
higher rate than Ri* (t) is equal to N — 1 — M(t). Hence, the 
total number of fractions used for channel probing within SDF 
algorithm at time t is given as follows: 

N p (t) = 1 + 1 + N - 1 - M(t) = N+1- M{t). (7) 

Note that N p (t) = N when all users are probed as in 
conventional Max Weight algorithm. 

We consider the following two functions: 



fs(Q(t))=E 



f m (Q(t))=E 



J2 (1 - Wp(t))W n (t)I n (t)\Q(t) 



n6Af p (() 



Y,(l ~ WWn(t)In(t)\Q(t) 



n£Af 



where the expectation is taken with respect to the randomness 
of channel variations and scheduling decisions. Given Q(t), 
both Max- Weight algorithm with full CSI and SDF schedules 
the same user which has the maximum weighted rate at every 
time slot. Hence, the value of W n (t)I n {t) is the same for 
both functions, and the only difference between f s (Q(t)) and 
fm(Q(t)) appears in the number of users probed. Next,we an- 
alyze the performance of SDF algorithm in terms of achievable 
rate region by using the theorem given in lETl . 

Theorem 5: lF2Tl If for some e > SDF algorithm 
guarantees 

fs(Q(t))>(l + e)f m (Q(t)) 

for all Q(t), then SDF can achieve 1 + e fraction of the rate 
region Af. 



Theorem 6: SDF algorithm can support (1 + e) fraction of 
the rate region A/, where 

_ /?(E[M(t)]-l]) 

e - — i^m — • (8) 

Proof: By using (0, f s (Q(t)) can be rewritten as follows: 



E 



5^(1 -f3(N + l- M{t)))W n {t)I n {t)\Q{t) 



f m {Q{t))+E 



^CSM(i)-£)W„(i)I n (i)|Q(i) 



(9) 



Now, we consider the value of f s (Q(t))/f m (Q(t)) such 
that, 



f s (Q(t))/f m (Q(t)) = 

f m (Q(t))+E \j2 neK(t) (l3M(t) - /3)W n (t)I n (t)\Q(t) 



f m {Q(t)) 

E[En(PM(t) ~ P)Wn(t)In(t)\Q(t)} 
fm(Q(t)) 



(10) 



Note that fh(Q(t)) represents the hypothetical capacity re- 
gion, Ah since j3 — 0. Hence, by comparing fh(Q(t)) and 
fm(Q(t))> we can find fmni which the maximum fraction 
that can be supported. 

Lemma 7: The maximum fraction e max is given by, 

f3N 

Proof: The proof is similar to the proof of Theorem |6] 
and is omitted. ■ 

Note that according to Theorem [6] the performance of SDF 
algorithm in terms of rate region depends on E[M(t)]. if 
E[A/(i)] > 1, then e > 0. In other words, we can increase 
the rate region. Next, we determine the the value of E[A/(i)] 
when channels are homogenous and heterogenous. 

V. Performance of SDF with Homogenous 
Channels 

We begin with a Lemma that shows how to determine the 
value of E[M (£)] when channels are homogenous. 

Lemma 8: When channels are homogenous, E[M(£)] is 
given as follows: 



where, 



E[M(t)] 



E 



Y, WM(t) - /3)W n (t)I n (t)\Q(t) 

n£Af p (t) 



Pi 



P2CY,P^ 



k=2 



Ps(J2p^ 

k=3 



vl 



(N-l) (15) 



iV-l 



J2 (fim - /3)E 



m=0 



Y W n (t)I n (t)\Q(t),M(t) = 

nGAf p (t) 



x Pr[M(t) = m] 
Note that 



(11) 



Proof: Let E[M(t)\R t ,(t) = r k ,n = i*] be the con- 
ditional expectation when the user which has the maximum 
queue length and its CSI are given. Note that for homogenous 
channels, the following equality holds, 

E[M(t)\Ri,(t) = r k ,n = i*] = E[M(t)\R,,(t) = r k ], (16) 

and E[M(i)] is determined as follows: 



Y W n (t)I n (t)\Q(t),M(t) = 

neN p (t) 

Hence, ( fTTT i can be rewritten as follows: 

N-l 



frn(Q(t)) 

1 - j3N 



E[M{t)\ 



53E[M(t)|Jfc.(t)=rfc]Pr[#,-.(t) 



r k 



fe=i 



m = ^^ J2 [Wm ~ /*)] Pr[M(t) = rn] 
Pf m (Q(t))(E[M(t)]-l}) 



1-/3 N 



Thus we have, 



fsW))/MQ(t))>i + WmW-il) (12) 

Hence, the proposed algorithm can support (1 + e) fraction of 

the rate region A/ where e = i-bn — • This com pl etes 
the proof. ■ 

Let us define 



When Ri*(t) = n, then SDF determines the user which has 
the maximum weighted rate by only using one fraction of 
time slot with probability pi, i.e., Pr[i?j*(£) = rj,] = pi. If 
this event occurs, the other N — 1 users do not report their 
CSIs with probability 1 and E[M(t)\Ri* (t) = n] = (N - 1). 
Similarly, when Ri*(t) = r<x, user j ^ i* does not report 
its CSI if Rj(t) < r 2 , and Pi[Rj(t) < r 2 ] = EJLjPk- 
Hence, with homogenous channels, E[M(t)\Ri* (t) = r 2 ] = 
(Sfc=2 Pk){N — 1). Now, for a given any value of L, we give 
a general formulation for E[M (£)] as follows, 



E[M(t)] 



Pi +P2(Yp^ +P^Y Pk ^ + 



vl 



k=2 



fe=3 



(N-l) 



f h (Q(t))=E 



J2 W n (t)I n (t)\Q(t) 

nGA/"(t) 



(13) 



Next, we investigate the case when the channels are ho- 
mogenous and uniformly distributed such that p k = j- for all 

k. 



Lemma 9: When channels are homogenous and uniformly 
distributed, E[A/(£)] is given as follows: 



E[M(t)] 



1 , L(L-l) 
L 



2L 2 



(N-l) 



1 1 

2 + 2L 



(iV-1) 

(17) 

Proof: If all channels are uniformly distributed, then pt = 
j-. In that case, by using (fl~5T > E[M(t)] is given by, 

E[M(t)} = 



fc=2 fc=3 



1 



Thus, we have, 

E[M(t)] = 



i + ^-l) + ^(X-2) 



1 

12 



(JV-1) 

(18) 
(N-l) 
(19) 



When we arrange ( fT9] l, we obtain 
1 L(L-l~ 



E[M(t)] 



L 



2L 2 



(N-l) 



1 1 

2 + 2L 



Theorem 10: When channels are homogenous and if N > 3 
we guarantee to expend the rate region, i.e., e > and the 
amount of increase in rate region is equal to e given in (|8). 

Proof: The proof of the theorem can be done in two parts. 
In the first part, we show that E[M(t)] is minimum when 
channels are uniform. In the second part of the proof, we 
show that if TV > 3, then e > when channels are uniform. 
We begin by proving the first part of the proof. 

We first show that E[A/(i)] is jointly convex function of 
P\ , P2 , ■ ■ ■ , Pl and then, the minimum of this convex function 
is achieved when channels are uniform. 

Lemma 11: E[A/(i)] is jointly convex function of 
pi,P2,---, PL- 
Proof: By using p 1 = 1 — ^ k=2 Pk, me Hessian of 
E[M (£)] in (fT~5T > can be given as follows, 

/ 2 1 1 1 1 
12 111 
112 11 
1112 1 
11112 



H = 



\ 1 1 1 1 1 



... IX 

... 1 
... 1 
... 1 
... 1 


... 2) 



(N-l) 



Now, we show that H is positive definite matrix. Let x = 
[xi X2 Xz ■ ■ ■ Xl] be any vector and x G R L_1 . If x T Hx > 
then, H is positive definite matrix and E[M (£)] is convex 
function of p 1 , p 2 ,..., pl El- 



x T Hx = 






\i=i 



•?7 



(JV-1) >0. 



(20) 



Hence, E[M (i)] is jointly convex function of pi,p2, ■ ■ ■ ,Pl- 

■ 
Lemma 12: E[M (t)] has the minimum value when channels 
are uniformly distributed. 

Proof: We already showed that E[A/(i)] is jointly convex 
function of Pi,P2, ■ ■ ■ ,Pl- Hence, by using p l = l-Y.k=2Pk 
we have the following L — 1 linear equations, 

SE[M(t)} 



SP2 

5E[M(t)] 

Sp3 

5E[M(t)} 
Sp 4 



5E[M(t)] 
Spl 



-1 + 2 P 2 + (P3+Pi + ---+P L )=0 (21) 
-1 + 2 P3 + (p2+P4 + ...+pl)=0 (22) 

-1 + 2p 4 + (P2 +P3 +P5 + ■ ■ ■ +PL) = 

(23) 
(24) 

-1 + 2p L + (P2+P3+P5 + ■ ■ -+PL-l) = 



(25) 
(26) 



(N — I) with matrix notation, 



( 2 1 1 1 1 

12 111 
112 11 
1112 1 
11112 



1 \ 

1 
1 

1 
1 



( P2 \ 



P3 
PA 
P5 



( l \ 



\PL J 



(27) 



\ 1 1 1 1 1 ••• 2 / 
Solving this linear system, we have, 

Pk~Pl = 0,Vfc,Z k ^ I 
Hence, 



\1/ 



Thus, 



Pk =pi,Vk,l k ^ I 



Pk = -T-,Vfc. 



(28) 
(29) 
(30) 



Thus, when the channel distributions are uniform, E[M (t)] has 
the minimum value. ■ 

We now prove the second part of the theorem. We show that 
when channels are uniform, E[M (£)] is a decreasing function 
of L. 

Lemma 13: E[A/(t)] is a decreasing function of L. 
Proof: From dTTb . 

"1 1 " 
2 + 2L 



E[M(t)} 



(N-l) 



(31) 



Taking the derivative of E[M (£)] with respect to L yields that, 



dE[M(t)} 

dZ 



-1 

2L2 



(JV-1) <0. 



(32) 



Thus, E[A/(i)] is a decreasing function of L. 



Now, it is easy to see that in ( fTTI i, taking L — >• oo yields 
that 



2.5 



x10 



lim E[M(t)] = lim 

L— ^oo L— >-oo 



1 1 

2 + 2L 



TV-1 
(#-!) = __ . (33) 



In the limiting case, when TV > 3, E[M(t)] > 1. In 
addition, according to Lemma [13] if L is finite, E[M(i)] is 
still greater than 1 whenever TV > 3 since E[M(i)] decreases 
as L increases. We can conclude that when all channels are 
uniformly distributed, and when TV > 3, then E[M(i)] > 1. 
As a result, we guarantee to expend the rate region, i.e., e > 0. 
In addition, according to Lemma [TT] and Lemma [T2l E[TW(t)] 
has its minimum value when channels are uniform. Therefore, 
for homogenous channels, when TV > 3, then e > and rate 
region is expended. This completes the proof. ■ 

VI. Performance of SDF with Non-homogenous 

CHANNEL 

When the channels are not identical, then p n k ^ p m k> where 
n y^ m and Vfc. In addition, for homogenous channels, the 
following equality holds, 



E[M(t)|-R i .(t)=r fc) n 



E[M(t)|i2i.(t) = r fc ], (34) 



However, this condition cannot be hold when channels are 
heterogenous. Hence, E[M(t)] cannot be determined in a 
similar way of homogenous channels. We have the following 
results for heterogenous channels. 

Theorem 14: When channels are non-homogenous and as 
L goes infinity, i.e., L — > oo, or TV is sufficiently large, then 
e > i.e., we guarantee to expand the rate region. 

Proof: the proof is provided in Appendix [A] ■ 

VII. Numerical Results 

In our simulations, we model a single cell CDMA downlink 
transmission utilizing high data rate (HDR) [56] ■ The base 
station serves 20 users and keeps a separate queue for each 
user. Time is slotted with length T s = 5 ms. Packets arrive 
at each slot according to Poisson distribution for each users 
with mean A n . The size of a packet is set to 128 bytes which 
corresponds to the size of an HDR packet. Each channel has 
5 possible states with rates as given in Table [IJ 

TABLE I: Possible Physical Rates 



Rates 


n 


T2 


»"3 


U 


»"5 


kb/s 


1843.2 


1228.8 


614.4 


307.2 


76.8 



A. Homogenous Channels 

First, we evaluate the performance of SDF algorithm when 
channels are homogenous and both uniform and non-uniform 
distributions. 
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Fig. 1: Performance of SDf algorithm with Homogenous and Uniform channels. 



1) Uniform Channels: In uniform case, p n = 0.2 for all 
users since there are 5 channel states and the channel state 
distributions are identical. Figure Q] depicts the average total 
queue sizes in terms of packets vs. the overall arrival rate 
when j3 = 0.01 and f3 = 0.02. The maximum supportable 
arrival rate is achieved in hypothetical case where the probing 
cost is zero, i.e., /3 = and the minimum supportable rate 
is achieved when full CSI is obtained and (3 = 0.02. When 
j3 = 0.01 and SDF algorithm is applied, we are closest to the 
hypothetical rate region. Clearly, SDF algorithm outperforms 
the full CSI Max- Weight algorithm in terms of rate region for 
both p = 0.01 and /? = 0.02. 

2) Non-Uniform Channels: Here, we investigate the perfor- 
mance of SDf algorithm when channels are identical but the 
channel state distributions are not uniform. In this case, the 
channel state distribution is given in Table HTl 

TABLE II: Channel state distribution 





Pi 


P2 


P3 


Pi 


P5 


probabilities 


0.3 


0.3 


0.2 


0.1 


0.1 



Figure [2] depicts the average total queue sizes in terms of 
packets vs. the overall arrival rate when channel state distribu- 
tions are uniform and as given in Table [TT] which is defined as 
non-uniform and p = 0.02. The maximum supportable arrival 
rate is achieved in hypothetical case the minimum supportable 
rate is achieved when full CSI is obtained. It is easy to see 
that the maximum supportable rate achieved by SDF algorithm 
when both uniform and non-uniform channels are the same. 

B. Heterogenous Channels 

Here, we investigate the performance of SDF algorithm 
when channels are neither homogenous nor uniform. For the 
heterogeneous case, we divide the users into four groups where 
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Fig. 2: Performance of SDf algorithm with Homogenous and Non-Uniform channels. 



Fig. 3: Performance of SDf algorithm with Heterogenous channels. 



TABLE III: Channel state distribution with heterogenous channels 





Pi 


V2 


P3 


P4 


P5 


Group 1 


0.3 


0.3 


0.2 


0.1 


0.1 


Group 2 


0.1 


0.1 


0.2 


0.3 


0.3 


Group 3 


0.25 


0.15 


0.1 


0.25 


0.15 


Group 4 


0.15 


0.25 


0.25 


0.1 


0.25 



x10 



there are 5 users in each gruop. The channel state distributions 
are given for each group in Table |TTT] 

Figure [3] depicts the average total queue sizes in terms of 
packets vs. the overall arrival rate when channel state distri- 
butions are given as in Table iLUl and j3 = 0.02. The maximum 
supportable arrival rate still is achieved in hypothetical case 
whereas the minimum supportable rate is achieved when full 
CSI is obtained. It is easy to see that the maximum supportable 
rate achieved by SDF algorithm when uniform, non-uniform 
and heterogenous channels are the same. 

1) Non-Stationary Channels: Next, we demonstrate the 
performance of SDF algorithm over time-correlated and non- 
stationary channels. The channel between the BS and each 
user is modeled as a correlated Rayleigh fading according 
to Jakes' model with different Doppler frequencies varying 
randomly between 5 Hz and 15 Hz. We set BW = 1.25 MHz 
and SNR = 10 dB. Let H n [t) denote CSI of user n at time 
slot t. H n {t) is a random process which does not have a 
stationary probability distribution, i.e, the mean of the channel 
gain changes over time. Let h n (t) represent the realization 
of H n (t) at time t, n G {1,2,..., N}. Then, the maximum 
number of bits of a user may transmit is given as, 

R n (t) = T s (l - (3N p (t))BW\og 2 (1 + SNR x h n (t)) (35) 




10 15 

mean overall arrival (packets/slot) 



Fig. 4: Performance of SDf algorithm over non-stationary channels. 



where BW is the bandwidth of a channel. Similar to the 
previous scenario, Figure |4] depicts the average total queue 
sizes in terms of packets vs. the overall arrival rate when the 
channels are non-stationary and f3 = 0.02. The maximum and 
the minimum arrival rates while keeping the queues stabile 
are achieved in hypothetical and full CSI cases, respectively. 
As the overall arrival rate exceeds 14 packets/slot queue sizes 
suddenly increase with full CSI Max-Weight and the network 
becomes unstable. However, SDF improves over full CSI Max- 
Weight by supporting the overall arrival rate of up to 16 
packets/slot. In other words, SDF algorithm is able to sustain 
14% more traffic than the full CSI Max-Weight algorithm. 
Thus, SDF appears to stabilize a larger range of arrival rates. 



VIII. Conclusion 

In this paper, we have considered the joint scheduling and 
channel probing problem in a single channel wireless downlink 
network. We have developed a low complex and dynamic 
feedback algorithm named SDF. We have shown that SDF 
algorithm can support 1 + e fraction of full achievable rate 
region. Then, we have proved the sufficient condition for 
e > 0. Our simulation results demonstrated a significant 
performance gain of the proposed algorithm compared to 
the case when the full CSI is obtained. In this extended 
abstract, the analytical and simulation results are given by 
assuming a single channel wireless network. In the full version 
of the paper, we will provide the result by considering the 
multichannel wireless system, i.e., OFDM networks. 
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Appendix A 
Proof of Theorem [T41 

A. SDF with Non-homogenous channel 

When the channels are not identical, i.e., p n k ^ p m k> where 
n / m and Vfc. We define the \n such that n = i*, i.e., user n 
is the user with maximum queue size at a time. In addition, we 
define tpf. such that Ri* (t) = r k , i.e., the user with maximum 
queue sizer is at channel state k at time t. Hence, 



Pr[x«] =Pr[n = i*], 
Pi[<Pk}=Pr[Ri*(t)=r k }. 

Then E[M(£)] can be found as follows: 

N L 

E[M(t)] = 53£)E[M(t)|xn, <Pk]Pr\Xn, <fih] 

n=\ fc=l 



(36) 



Hence, 








E[M(t)} = 






( 


N L 


N \ 


Pr[xi] pn(iV- 


- l)+Pl2C^2^2Pnk) + ■ 


• + PilC^Pul) j 


V 


n=1k=1 


n=2 / 


/ 


N L 


\ 


+ Pr[x2] P2i(JV- 


-l)+P22(^^Pnfc) + . 


• +P2LC^2PnL) 


V 


n=l fc=2 
n^2 


n=l 

n#2 / 


/ 


N L 


N \ 


+Pr{xs] 


P3i(N- 


-i)+p3 2 (J2J2p^) + - 


■ +P3LC^2PnL) 




\ 


n=l k=2 
n/3 


n=l , 
n#3 / 



N L 



N 



■Pr[xjv] \pNi{N-l)+p N 2(J2J2 Pn ^ 



Pnl 



(^ PnL) 



ra=l k=2 
n=jtN 



n=l 



Note that Pr[x„] > p" q lin for all n, where < p" q lin < 1. 
Hence, a lower bound on E[M(i)] can be given as follow, 

E[M{t)] > Np™ in [p min (N - l)+p rnin (p min {L - 1)(7V - 1) 
+ p min (p mm {L - 2)(/V - 1) + . . . 

+ p min (p mm (N -1)] 

By rearranging, we have, 



E[M(i)] > Np 



in in 

q 



(n - 1) P mm + ( P mm y 



L(L-l) 



10 



Therefore, we can guarantee to achieve a larger capacity region 
when the following condition is satisfied: 

N(N- 1) > 



n 



p" 



(yjmin\ 



L(L-l) 
2 



Clearly, L — > oo or TV is large enough, this condition holds 
and e > 0. This completes the proof. 



